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Abstract 

We study the radiative leptonic —)• ytn decays in the nonrelativistic QCD effective field 
theory, and consider a fast-moving photon. As a result the interactions with the heavy quarks 
can be integrated out, and thus we arrive at a factorization formula for the decay amplitude. We 
calculate not only the relevant short-distance coefficients at leading order and next-to-leading order 
in Us, but also the nonrelativistic corrections at the order |vp in our analysis. We find out that 
the QCD corrections can sizably decrease the branching ratio and thus is of great importance in 
extracting the long-distance operator matrix elements of B^- For the phenomenological application, 
we present our results for the photon energy, lepton energy and lepton-neutrino invariant mass 
distribution. 
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I. INTRODUCTION 


The search for new degrees of freedom can proceed under two distinctive directions. At 
the high energy frontier, new particles have different signatures with the standard model 
(SM) particles, and measurements of their production may provide definitive evidence on 
their existence. On the other hand, it is likely that low energy processes will be influenced 
through loop effects. Rare decays of heavy mesons, with tiny decay rates in the SM, are 
sensitive to the new degrees of freedom and thus can be exploited as indirect searches of 
these unknown effects, for a recent review see Ref. [1]. 

The Be meson is the unique pseudo-scalar meson that is long lived and composed of 
two different heavy flavors. Since this hadron is stable against strong interactions, its weak 
decays provide a rich phenomena for the study of CKM matrix elements, and also a platform 
to study the effects of weak interactions in a heavy quarkonium system HE]. In the 
past decades it has received growing attentions since the first observation by the CDF 
collaboration [1]. This can be particularly witnessed by the recent LHCb measurements of 
the Be lifetime [HIE], the decay widths of Be — )■ J/'07r and Be —)■ J/'tpi.i' [3 E], and various 
other decay modes [9HI2]. One may expect that more decay channels of Be can be measured 
by the LHCb, ATLAS and CMS experiments [T3lfT5] . 

On theoretical side, various approaches have been applied to calculate the decay width 
of Be decays [T6H52|, but most of them are phenomenological. Since both constituents of 
the Be are heavy and can only be treated nonrelativistically, an effective held theory can be 
established [53]. Taking the Be —)■ J/iIjIu as the example, one may derive the conjectured 
non-relativistic QCD (NRQCD) factorization formula for its decay amplitude: 

A{Be ^ J/^) oc C,,{Q\Of\Be) X (J/^|Of |0), (1) 

f f 

where the (9/j are constructed by low energy operators. The short-distance, or hard, con¬ 
tributions at the length scale are encapsulated into the coefficients Cij that can be 

computed in perturbation theory. 

The long-distance, or soft part of, matrix elements have to be extracted in a nonpertur- 
bative approach, for instance the Lattice QCD simulation, or constrained by much simpler 
processes for instance the annihilation modes Be —)■ iit and Be —)■ yf'P. However, the useful¬ 
ness of the Be —)■ £i^ is challenged by two aspects. Firstly its decay rate is given by 


T{Be ^ iue) 


^\ycb\"fkrn[ 


.3 
Be 


m 


Be 



( 2 ) 


in which the suppression factor arises from the helicity flip. As a result, the Be —)■ 

and Be —)■ ePg have tiny branching fractions that may be out of the detector capability at 
the current experimental facilities. Secondly, there is only one physical observable, namely 
the decay rate, and thus the Be —>■ f'P is not capable to uniquely determine all, typically 
more than one when relativistic corrections are taken into account, long-distance matrix 
elements (LDMEs). 

On the contrary, the Be —)■ 7^P can provide a wealth of information [53HHH], in terms of 
a number of observables ranging from the decay probabilities, polarizations to an angular 
analysis. It is interesting to notice that the counterpart in B sector, B —)■ y^P, has been 
widely discussed towards the understanding of the B meson light-cone distribution ampli¬ 
tudes [591 - I63] . The small branching fraction of Be —)■ y^P can be compensated by the high 
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FIG. 1: Leading order Feynman diagrams for the radiative leptonic Be —)• decay in the SM. 

The lepton /i can also be e or r. The photon emission from a virtual VF-boson shown in the second 
panel is suppressed by 1/m^ compared to the other contributions. 


luminosity at the ongoing hadron colliders and the under-design experimental facilities. The 
main purpose of this paper is to explore the t 7 ^z/ at next-to-leading order (NLO) in 
as and in |vp, which shall catch up the progress in the Be —)■ [SSIEI]. For the leptonic 


decay constant, the two-loop calculation is also available in Ref. [65] . 

The rest of this paper is organized as follows. In Sec. |T| we will derive the formulas 
for various partial decay widths of Be —)■ yt'P. Sec. uni is extensively devoted to the next- 
to-leading order calculation. We will discuss the phenomenological results in Sec. [IV} We 
summarize our hndings and conclude in Sec. |V| We relegate the calculation details to the 
Appendix. 

In the SM, leading order (LO) Feynman diagrams for the Be —)■ yt'P decay are shown in 
Fig[^ The photon emission from a virtual LF-boson is suppressed by 1/m^ compared to 
other contributions, and thus the second diagram in Fig. [T] can be neglected. Integrating 
out the off-shell LF-boson, we arrive at the effective electro-weak Hamiltonian 


G 

HeE = -^VebC-i^il - 75 ) 6 / 7'"(1 - 75 )// h.c., 


(3) 


where Veb is the CKM matrix element. The decay amplitude, matrix element of the abo 


II. Be -llv 

ve Hamiltonian between the Be and yf'P state, 

A = (7rz/|iLefr|5e) (4) 

is responsible for the process Be —)■ ^Iv. 


A. Differential decay widths 

Since there is no strong interaction connection between the leptonic and hadronic part, 
the decay amplitude can be decomposed into two individual sectors: 


A= (0|c7^(l -75)&|5c) X (7/ - 75)i/|0) 


-K(7|c7m(1 - l5)b\Be) X - 75)u, 


(5) 
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with the matrix elements encoding the hadronic effects: 


(0|c7^(l - 75)&|5 c), (7|c7/.(l - lb)h\Bc). 

The first one dehnes the Be decay constant 

(0|c7^75&|5c(pbJ) = 


(6) 

(7) 


while the Be ^ 'y transition is parametrized by two form factors: 


(7(e,/c)|c7^6|5c(pBj) 
( 7 (e, A:)|c 7 ^ 756 | 5 c(pbJ) 


-e^ 
PBe ■ k 

ieA{L‘^) i 


^{lUpCT^ FBc ’ 

e* -k 

^ ^PB.-k 


ie 

Pb, ■ k 


IBcPBcIiPBc 


( 8 ) 

( 9 ) 


with the momentum transfer L = — k. Here and throughout this work we adopt the 

convention rjpg above equations are similar with the parameterization of the 

H —)■ 7 form factors as given in Ref. [HO]- The last term in Eq. (|^ that is proportional to 
the Be decay constant has been added in order to maintain the gauge invariance of the full 
amplitude [HTJEH], and see appendix [A| for a derivation. 

Substituting Eqs. @ into Eq. ([^, we obtain 


A 



VebefBe^iYil 




[1 + a(s/)] 



PBe-k ) 


iv{si) 
Pb, ■ k 








( 10 ) 


where si = and terms due to lepton mass corrections have been neglected. Apparently, 
this expression is gauge invariant. For the sake of simplicity, we have dehned two abbrevia¬ 
tions in the above ^ 


^ A{si) R(sO 


aiB) ^ 


fBe 


v{si) = 


Jb, 


( 11 ) 


In terms of the decay constant and form factors, the differential decay width for the 
Be —1is given as 




dEkdEi 




QArriB^T^^ 

4:7l‘^xl 


(1 - Xk) X 


{xl + 2xk{xi - 1) + 2{xi - 1)^) 


+2a {{v + l)xl 2{v l)xk{xi - 1) + 2{xi - 1)^) 2vxk{xk + 2xi - 2) 
-f-n^ [xl + 2xk{xi - 1) 2{xi - 1)^) + xl + 2xkXi - 2xk + 2x1 -4:Xi + 2 


(12) 


One shall distinguish the form factor v from the relative velocity v to be defined in the following. 
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where Xk = 2Ek/mB^ and y = 2Ei/mB^, and Ek and Ei is the energy of the photon and 
charged lepton in the Bf. rest frame, respectively. One can integrate ont the Ei and obtain 


dr _ aemfBj^cbl Gj,m%^Xk{l - Xfc)((l + g)^ + v^) 
dEk 12 tt‘^ 

The differential distribntions can also be converted to 


(13) 


£r 


dsid cos 0i 


|KbPaem/s,G|(l - Xk)^ 


Xt 


{xl + 2xk{xi - 1) + 2{xi - 1)^) 


ms, - Si 

32mB,n‘^ 

+2a [{v + l)xl + 2{v + l)xk{xi - 1 ) + 2{xi - 1 )^) + 2vxk{xk + 2xi - 2 ) 
[xl + 2xk{xi - 1) + 2(xi - 1)^) + xl + 2xkXi - 2xk + 2x1 -4x^ + 2 


(14) 


nsing the relation: 


Ek 


El 


- Si 


2mB, 

[i^k + - i^k 

Be 


Si) cos 0i] . 


(15) 

(16) 


The 6i is the polar angle between the lepton i flight direction and the opposite direction of 
the Be meson in the rest frame of the iui pair. Likewise one can integrate ont the 6i 


dT _ Q:em/|jV'c6pG|,(m|^ - si)si{{l + af + 
dsi 247r‘^m^^ 


B. NRQCD factorization 


The factorization properties for the B^ —)■ yf'z/ depend on the kinematics of the photon. 
In this work, we will not stndy the soft-photon contribntion as disenssed in B decays [5U] . 
and leave it for fnture work. In the region where the photon is a collinear (fast-moving) 
object, its interaction with heavy qnarks is highly virtnal and thns should be encoded in the 
short distance coefficients. In the NRQCD scheme, we only need retain those color-singlet 
operator matrix elements that connect the Be state to the vacuum. To the desired order, 
one expects the following factorization formula: 


/n, 

V 

A 



Co(Oixl^bjBe(p)} + ^(Ojxl i^bjBe(p)} + O(v^) , 

— (Ojxl^biBe(p)} + ^(Olx); ^biBe(p)} + O(v^) 

tub^ \ z j 

— (0L\oV'i|Sc(p)> + -T( 0 L\i f-^Tj i^i|B.;(p)) + ©(v*) 

jub^ rfig \ Z / 


(18) 

(19) 

( 20 ) 
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where v denotes half relative velocity between the charm and bottom quarks in the meson, 
and C 2 ^’^ are the dimensionless short-distance coefficients that can be expanded in 
terms of the strong coupling constant We shall calculate the one-loop corrections to 
the but give only the LO results for since the latter ones are already power- 

suppressed. 'ijjQ and Xq represent Pauli spinor helds that annihilate the heavy quark Q and 
anti-quark Q, respectively. Besides, one need note that the state \H{p)) in QCD has the 
standard normalization: {H{p')\H(p)) = 2Ep{2n)^S^{p — p'), while an additional factor 2Ep 
is abandoned in the nonrelativistic normalization where (i?(p')|iP(p)) = (27r)^(5^(p — p'). 


III. NEXT-TO-LEADING ORDER CALCULATION 
A. Kinematics 

Let Pi and p 2 represent the momenta for the heavy quark Q and anti-quark Q\ Without 
loss of generality, one may adopt the decomposition: 

pi = aPB,-q, ( 21 ) 

P 2 = Pb, + q, (22) 

where Pb^ is the total momentum of the quark pair, g is a half of the relative momentum 
between the quark pair with Pb^ ■ q = 0. a and (3 are the energy fraction for Q and Q' in 
the meson, respectively. The explicit expressions for all the momentum in the rest frame of 
the Bf. meson are given by 


n. 

— {El -|- id2, 0), 

(23) 

q^ 

= (0,q), 

(24) 

Pi 

= (^1,-q), 

(25) 

P 2 

= (i^ 2 ,q). 

(26) 


In the rest frame, the meson momentum becomes purely timelike while the relative momen¬ 
tum is spacelike. One can obtain the relations a = — q^) and 

/d = 1 — a with the on-shell conditions Ei = \/ml — E 2 = \/ml — g^, and g^ = —q^. 


B. Convariant projection method 


In the following calculation, we will adopt the covariant spin-projector method, which 
can be applied to all orders in v. 

The Dirac spinors for the system may be written as 


Ub{pi,\) 


El mb 
2Ei 


6 

Ei+rrif, 


6 


(27) 


^ Throughout this paper, we shall use the superscripts (0) and (1) to indicate the LO and NLO contributions 
in as and the subscripts 0 and 2 to denote the LO and NLO contributions in the velocity. 
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Vc{p2, A) 


(28) 


E 2 + rric 


2Eo 


E2+mc^^ 

6 


? 


where is the two-component Pauli spinors and A is the polarization parameters. It is 
straightforward to derive the covariant form of the spin-singlet combinations of spinor bilin- 
ears: 


no(g) = V Ai)uc(P2,A2)(^Ai^A2|00) ® ^1= 

^ ^ yEc 


Ai,A2 

i 


A\/‘2EiE20J 


/ / / , N i^Bc -\- El -\- E 2 , „ , , . Ic 

[oipB, -i + mb) - -^^^^2 -^ - mj ® , 


(29) 


with the auxiliary parameter u = \/Ei + mb\/E 2 + rric. Here Ic is the unit matrix in the 
fundamental representation of the color SU(3) group. 


C. Perturbative matching 


Due to the simplicity of the final state, one can directly match the QCD currents onto the 
NRQCD ones. To determine the values of cq and C 2 , we follow the spirit that those short- 
distance coefficients are insensitive to the long-distance hadronic dynamics. As a convenient 
choice, one can replace the physical B~ meson by a free cb pair of the quantum number 
so that both the full amplitude, ^[c6(^S'o^') —)■ and the NRQCD operator matrix 

elements can be directly accessed in perturbation theory. The short-distance coefficients 
Ci can then be solved by equating the QCD amplitude A and the corresponding NRQCD 
amplitude, order by order in For this purpose, we introduce a decay constant and two 
form factors at the free quark level: 


{0|c7^75&|c&(^4^^)) = 


(7(e,fc)|c7^6|c6(^4^^)) = 

ft ■ PBa 


( 7 (e, k)\c^^^A>\ch{^S\i"')) = ieA ( e* - M - ie — ^—tI^Pb^i^Pb, ■ e*. 

V PB,-k J PB,-k 

Analogous to ( I^19[20 ), one can write down the matching formula: 


^ = 4{^\xl'4’b\cbCsi^^)) + 

1 


‘-2 


AlcbCsll^)), 


V = 


A = 


rUb + rric 
1 

mb -|- me 


mb + ruc)^ 
Coi^lxlHcbCsj^^)) + 7 — 


-.v 


(mb + mc)^ 


Coi^lxlHcbCsfl^)) + 7-X- 

[mb + me) \ ^ J 


(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


where we have adopted the nonrelativistic normalization. 
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One can organize the full amplitudes defined in Eqs. ( ^31[32 ) in powers of the relative 
momentum between c and 6, denoted by q. To the desired accuracy, one can truncate the 
series at C>(q^), with the hrst two Taylor coefficients. We will compute both amplitudes 
at LO in in subsection HID, and the calculation at NLO in will be conducted in 
subsection IIIIEl 

The NRQCD matrix elements encountered in the above equations are particularly simple 
at LO in a*: 




(36) 


where the factor \/2Nc is due to the spin and color factors of the normalized c6(^S'q^^) state. 


The computation of these matrix elements to 0(as) will be addressed in subsection IIIF 


D. Tree-level amplitude 

Adopting the above notation, one can easily obtain the tree-level amplitude for the decay 
constant 


( 0 |c 7 /, 756 |c 6 (^ 4 ^^))(°) = Tr [no(g)7M75] 

(El + mb) (Ea + 


= 


2\/ EiE2(Ei + mfe)(E2 + rUc) (Ei + Ea) 








where the terms have been omitted and 


fBred — 


rribmc 
mb + me' 


is dehned as the reduced mass of the cb system. 
The vector current is similarly evaluated as: 


(37) 


(38) 


/ I- c'[d^\(o) T" rrr z' ^ ~ 1^2 + "ic) „ r„ . ^ i{'^i — ^ + mb) . 

(7|c7^6|c6( ^ = TT[Uo{q)teec^ - -+ Tr[no(g)7M7-Tw- \ 


{k - pa)^ - ml 

e\/2Nc , Cc 


(pi — /c)2 — ml 


+ 


eb 


dtci/EiEa Ea/c ■ Pb^ + Ek ■ q Eik ■ - Ek ■ q 

X IEbc^^upa^*’'k^p^^ + E(Ei + Ea + mb — mc)^ij,vp(jk*'^k^q° 


(39) 


We have introduced the abbreviation E = Ei + Ea, and Ebc = (Ei + mb){E 2 + me) + q^- 
Here Cc = 2/3 and eb = —1/3 is the electric charge of the c and b quark, respectively. 

One can perform the Taylor expansion of the amplitudes in powers of q^: 


A{q) 


^(0) + 


a^(o) 

dq^^ 


q=o q^ + 


1 aM(o) 


2! dqf^dq'" 


,=o q^q'" + ■■■■ 


(40) 

























Those terms linear in q should be dropped since this auxiliary momentum introduced at the 
quark level has no correspondence at the hadron level. In this paper, the (9(|qp) contribu¬ 
tions will be retained. In order to simplify the calculation in the covariant derivation, one 
shall use the following replacement: 




-)■ 


|qP 


D - 1 


+ 


Pb Pr 


)• 


(41) 


The result for the axial-vector current is a bit lengthy: 


(7|c7/.75&|c6(^4^^))^°^ 


* 




ec 

E2k - Pb, + Ek ■ 

q 



k-PB, 

Ebc + k ■ qE{Ei - E2 + mb 

- rUc) 


» 

Elk ■pB,-Ek-q 




2{Ei- 

- E 2 +mb- mc){E2e* ■ 

Pb, 

+ Ee* ■ 

■q) 


E2k ■ Pb, + Ek ■ 

9 



2{Ei- 

- E 2 +mb- mc){Eie* ■ 

Pb, 

- Ee* 

■q) 


Elk ■ Pb^ — 

2Ebc{E2e* ■ Pb, + Ee* ■ q) 

+PB,f.e, + Ek-q) 

2(EiEbce* ■ Pb, + Ee* ■ q{Ebc + g^) 
E{E,k ■ PB, - Ek ■ q) 


-PB,bLeb -ip/jp , - 

E[Eik-pB, 

EbcP • Pb, + Ee* -qjEi- E2 + mb 

E2k ■ Pb, + Ek ■ q 

Ebce* ■ Pb, + Ee* ■ q{Ei - E2 + mb 

+kfj.eb -- — - 

Elk ■ Pb, - Ek ■ q 


-kf,e, 


(42) 


In order to extract the A form factor, we only need to keep the term which corresponds 
to Feynman gauge e ■ Pb, = but we have explicitly checked the gauge invariance up to v^ 
order. 

The tree-level NRQCD matrix elements for the c6 have been given in Eq. (36), and thus 
the above results in Eqs. (^39[42) lead to the tree-level Wilson coefficients 


/.o 

Co = 

1, 

Jfl _ 


C 2 

8^2’ 

wo 

6c 65 

Cq — 

1 

to 1 
0 1 

1 

C 2 — 

-2 f ec(3z‘ 

Afl 

\ 

65 6c 

Cq — 

2 2z’ 


ee(3z2 + 2z + 11) ^ Cbillz^ + 2z + 3) 


48^3 


482;^ 


(43) 

(44) 

(45) 

(46) 

(47) 
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FIG. 2: Typical NLO Feynman diagrams for the radiative leptonic —)• decay in the 

SM. The other four diagrams can be easily obtained by interchanging the bottom and anti-charm 
quarks lines. 



~2 f ec[(3;2^ + 2^ + 11 ) + 82;(1 - z)mb/Ek] 

\ 48^;^ 

eb[(llz^ + 2^; + 3) - 8z(l - z)mb/Ek\\ 
48^2 ]■ 


(48) 


In the above results, we have dehned 2 ; = mc/rrib and z = 1 + z. cf’° means the LO of 
Wilson coefficient cf. It is interesting to notice that the Wilson coefficients depends 
on the energy of the emitted photon, which will induce nontrivial behaviors as will be 
demonstrated later. 


E. NLO amplitudes in QCD 


Typical one-loop diagrams for the QCD corrections to the Be —>■ decay are shown 

in Fig. In calculating the one-loop amplitudes, we use the dimensional regularization to 
regulate the ultraviolet (UV) and infrared (IR) divergence. 

The diagram (a) in Fig. [^contributes to the NLO decay constant: 


with 


Vl.a = \/^ 


■CpCis 

4:71 


1 2 6 In 2 ; 

V-1“ 2 -h 3 In —^ — 2 2ti -- 

^uv ^IR z +1 


U = 


V = 


1 

1 9 

V 

_ 1 

1 


T - In 

e/ij h 


q 


‘2,7Tlred 


(49) 


(50) 


We have introduced the abbreviation 

1 


^UV,IR ^UV,IR 


'jE + In dvr. 


(51) 
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The heavy quark field renormalization and mass term are given as 


1 - 


= 


Cpas 
dvr 
SmCrai 
dvr 


v^ + A + 31n-^ + d 
euv ^iR 

■ 1 , /x' d 

-h hi —- + — 


. O • (52) 

_euv ^ 3_ 

For the vector current form factor, the sub-diagram in Fig. gives out the corresponding 
contribution 

y^cetCpas. 1 d 5 -2 , „.2 1 0X..2 


V, = 


Annib 

22/2 


[-T, 


1 ,. ,2 (9=+ 1 )), 

^-ry + — -O2 i--ry;-C»3 

^iR y ~ z y z — z-^ z z [y^ — z^) 


y 


Z (2/2 — p 


’hi - —C4 - (1 - Z)C^ - (^2 _ y^)di], 

z 


V. = 


\/2NcebCFQ-s r 1 1 y'^ + z‘^ + Az + 3 z + y 


+ 


Anmb 2 euv 
(2^2 -|- 3 ^ _ ^ _ |/2(2z - 1 - 3) 

2z ( 2/2 — P) 


z2 _ y2 


+ 


2;2 - 1/2 


& 1 - 


^ (32/2 - ;g2 ^ 1 ) 

2^ (52 — 1/2) 


^3 + (^ + 2 /^ -2;^)c4], 


Vrf = 


y/P^ebCpas . 1 1 y'^ - z‘^ + Az + 3 y^ - z"^ + z + 2 

-'-^-^-^--3-O 2 


dvrmb 2 euv 

y'^ — z"^ + Az + 5, , 

03 + ClJ, 


z2 _ y2 


z2 _ y2 


2 ( 2/2 - z‘^) 

_ ^flNcebCpOis^-y"^ + z^ + 8z + 7 ^ z 

V — ; I 7 rT~Zn fT\ “r 


'f _ Z ( 2 /^ - ^2 + 1 ) 

dvrmb '■ 2(52 — 7 / 2 ) ' 2 ( 2/2 — z5) 2 {y‘^ — zz) {y‘^ — P 


2 _ 2 / 2 ) 

(2;2 + Qz + 1) P + y‘^ — 22/2 (^2 -| 4^; + 3) 


2 ( 2 /' 


( 2 /^ 






(53) 


where the auxiliary functions hi, Ci, and di are defined in Appendix [Bj 

The counter-mass terms and wave function renormalization corrections give: 


^CT-m — 


^CT-F — 


y/2NcebCFas ^ 3z 
Anrub 

\/2NcebCFas 


1 , 2 x 2 q 

L 2 ~2 V --^ -2 + f) ’ 

2/2 - 2;2 3 

, 1 11 3, /x2 

+ X- 1 " In- 2 ^ • 

2 zmi 


Anmb 'em ' 2 euv ' 2"^ zm^ ' '' 

For the axial-vector current form factor, the sub-diagram has gauge-dependent contribu¬ 
tions, however, the summed result is gauge-invariant. We will show the detail in Appendix [C| 


F. NLO amplitudes in NRQCD 


The NRQCD Lagrangian can be derived by integrating out the degrees of freedom of 
order heavy quark mass [3H] : 

/ Y)2 \ j^4 ^ 

>^^NRQCD = ^ j ^ 

+ - 9s^ ■ ■ (D X fif^E - QsE X D)i/> 

-|- (-0 —)■ iPx* 1 — A^) -|- Aight • (55) 
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The replacement in the last line implies that the corresponding heavy anti-quark bilinear 
sector can be obtained through the charge conjugation transformation. £iight represents the 
Lagrangian for the light quarks and gluons. The coefficients c^, cp, and cs have perturbative 
expansions in powers of which can be written as c* = 1 -|- 0{as). 

The matrix element of the cb to vacuum at NLO can be written as 

= y^^ 2n2\v\ 

This is in agreement with the results in Ref. uni. 


eiR /i^ 


( 56 ) 


G. Determination of cp Matching QCD to NRQCD 

Up to tts and v^, one can expand the decay constant and form factors as 

u = ci’“{0|xJVi|fK‘4")>‘“’ + d'CJIxJii’iIfK'si"))'"’ + d'CJIxSV'iIcK'si"))*'’ 

Jfi 




V = 


1 




mb + me 


^2 




A = 


mb + me 




+4’"{o|\'J4°’|cK‘4‘')>“> + 




(m. + m.)4'>l4(-2^)4.I^K'4‘')>'“']. 


( 57 ) 


( 58 ) 


( 59 ) 


Matching the QCD results onto the NRQCD, one can obtain the UV and IR finite short- 
distance coefficient 


/4 

Cn’ = 


3C,a. 


dvr 


1 + z 


( 60 ) 


VI CpOls r r, M 
Cn = — —{e^ln — 


An 


^ {—?>zzz + 2y^) + z'^ + y'^{?>z — 1) ^ y‘^ — 2zz , 


2z 


zml 2 {y^ — zz) {y^ — z^) 
2 


4 [z'^ — y‘^z) 


2z — zz) 


_ - ^_ --‘Wh + -zz‘^ + z‘^ + ^y‘^z-y‘^ 

A^y‘^ — zz z — y z + y z ^ 2z{y‘^ — P) ^ 


—z — y'^z + z^ + z^ 
zz 


y'^z - z^ + 2 z + I . ,, . 2 -25 7 1 

Cl H-C2 + {z — l)c3 + {y — z )di] 


+{eb —t ——)■ -, 1 / —)■ -)}, 
z z z 


( 61 ) 


12 



























s,(GeV 2 ) 


i,(GeV^) 


FIG. 3: Dependence of short-distance coefficients on the si. The solid line denotes the 

coefficient the dotted line is the coefficient from relativistic corrections, and the 

thick curve is the coefficient from as corrections. 


A,1 

Cn = 


CpCHs 

4:71 


Gb 


lnA.+ 


zm^ 2 (|/2 - zz) (|/2 - py 


+ {z{z{3z + 23) + 5) + 1)5=^ + /) + 
-3)z^ + y\3z-l))+ 

+ {3{z-l)z-2)z^) - 


4z (r/^ — PY 


;{y‘^{z + 11)5 — y‘^{z{5z -h 34) + 5)z‘^ 
{-2y‘^{z - 3)P - + Uz 


4z — zz) — z"^) 
bs 

2z{zz — y'^){y‘^ — PY 


(r/2(2/2(3^ + 7)_(2^ + 3)(3^_l)5) 

[y^ (13^^ - 2;^ + l) 5^ - 2(3;^^ -h z)P 


- 8 .^ - 6 . + 2 )) - + + + + 

^ 2z{7P-Pf 

y‘^ {y‘^{—z) + z‘^{2z -|- 5) — 3) — (z — l){z{z -|- 4) — 1)5^ 


zz {P — y"^) 


Cl 


((z - 2)z(z + 4) + 1)5" + y^ (z (y^ - 2z(z + 2) + 3 ) + 3) 
z(y-z)(z + y) 

(Z — 1) (— -h Z^ — 1) , 2 2 

+ ^- -C^ + (V + + 4z - l)c?i 


-C 2 


y2 _ ^2 

-\Cb —)■ —,z —)■ -,y -A- -^ 
z z z 


(62) 


Note that the scale dependent term in the brace of Eqs. (61) and (62) will be cancelled each 


other, the residual dependence only lies in the strong coupling constant. 


IV. PHENOMENOLNICAL RESULTS 

The input parameters are adopted as [7T]: = 6.2756GeV; Gp = 1.16637 x 

10“^GeV“^; a = 1/128; for the GKM parameters, we adopt \Vcb\ = 0.041. For the heavy 
quark mass, we adopt rrih = 4.8GeV and rric = l.SGeV [IH]. The i?c"Kieson lifetime is using 
the latest measurement by the LHGb Gollaboration, i.e. = O.SOps ISIIB]- 
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We first present numerical results for the decay constant fs^- 


./>0 _ 

X 

-3.8, 

•2 

8x2 

./-I 

3C pOig 

(2± 

0 ~ 

An 


1-z 



—0.44 X a<j. 


The strong coupling constant at the Z-boson peak is mi 


a.(mz) = 0.1185 ±0.0006, 


which corresponds to 


otsijnb) = 0.218, asirric) = 0.368. 


(63) 


(64) 


(65) 


With these values, one can see the corrections can reduce the decay constant by approx¬ 
imately 9.5% — 16.2%. 

To estimate the size of (!2(|vp) effects, one requests the size of non-perturbative LDMEs, 
for which we use Buchmiiller-Tye (B-T) potential model [72] : 

(0|xi'!/’b|5c(p)) 

V’b|5c(p)) 

For an estimate of q^, one may make use of the relative velocity. Using the heavy quarks 
kinetic and potential energy approximation [53], we have 

|v| ~ as{2mred\y\) ■ (68) 


= \ — I - 0.884GeV^/2 , (66) 

V 27r 


- q^{o|\'SV’i|Bc(p)). 


(67) 


Choosing mb = 4.8 GeV and rUc = 1.5 GeV, and using two-loop strong coupling constant, 
we get 

|v|^/^ 0.267, 0.108, |v|^^ 0.186. (69) 

For a value ~ 0.186, we have 

~ 0.9718GeV^ (70) 


As a result, the decay constant will be further reduced by about 

For the short-distance coefficients for Be ^ j transition form factors V and A, our results 
are shown in Fig. ^ The solid line denotes the leading-order coefficient the dotted 

line correspond to the coefficient from relativistic corrections, and the thick curve is 

the coefficient from corrections. From these hgures, one can see the relativistic 

corrections give constructive contributions, but the 0{as) QGD corrections are destructive 
and thus have important consequences. Note that the factorization in Fqs. ( I^19[20 ) is 
valid only for a hard photon, while the soft-photon contribution needs special treatment [69] . 
Thus a cut-off on the photon energy should be introduced, however we have checked that 
the cut-off will not affect the results signihcantly in Tabs. [T| and |TTj 
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FIG. 4: The dependence of the branching ratio B{Bc —)• on the photon and lepton energy. 

The dotted line denotes the leading-order result, the dashed line is the result with relativistic 
corrections, the blue line is the result with QCD corrections, and the thick curve denotes the total 
results with both the QCD and relativistic corrections. 



FIG. 5: Similar with Fig. but for the si dependence. 


With the estimated long-distance matrix elements, results for differential distributions are 
given in Figs. |^and[^ where the QCD and relativistic corrections are shown respectively. 
The integrated branching ratios of —)■ 'jii^ and —)■ ii^ are presented in Tabs. and 
[nj Ignoring the lepton mass, the branching ratio of Be —)■ yez/g is identical to that of 
Be —)■ The LO results are in agreement with Ref. [HTl - IHS] with the same input 

parameters. From the calculation, one can see that both the QCD and relativistic corrections 
give destructive contributions to the process Be —)■ iiy. However, relativistic corrections 
produce a constructive contribution to the Be —)■ Our results have demonstrated that 

the QCD and relativistic corrections are mandatory towards a more accurate extraction of 
the value of LDMEs for Be system. 
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TABLE I: Branching ratios of Be —)• 'yf.v and Be —>■ d-v. Here tb^ = 0.50ps, and we vary the heavy 
quark masses with mb = 4.8 ± O.lGeV and me = 1.5 =F O.lGeV. 


Channels 


Tree-level 


|v|^-corrections QCD corrections This work (NLO) 


Be- 

->• TU-r 

2.90 X 10-2 

-0.54 X 10-2 

—0 

^•^*^-0.04 

X 10-2 

1 

-‘-•^^-0.04 

Be- 


12.10 X 10-5 

-2.25 X 10-5 

9 00+0.14 

X 10-5 

7 cq-rO.M 

Be- 

eue 

2.82 X 10-9 

-0.53 X 10-9 

-0.54l9;01 

X 10-9 

1 yc+U.U^ 

' ^-0.04 

Be- 


10 . 49 + 2-27 X 10-5 

5.46^+^^ X 10-5 


X 10-5 

o qo+1.65 

o.zo i 33 


X 10“2 
X lO-*^ 
X 10-9 
X lO"*^ 


TABLE 11: Branching ratios of Be —)• and Be —t- ii' compared with other theories or models, 

including Lattice QCD (LQCD), Light front model (LEM), Constituent quark model (CQM). Here 
tb^ = 0.50ps is adopted. 


This work LQCD [71] LEM |56| CQM [58| Ref. [75] Ref. [55] 


10^B{Be 

TUr) 

1 or\+0.03 
-*-■^^-0.04 

2.12 

1.52 

1.44 

1.8 

1.6 

W^B{Be 


7 cq-ro.u 
' ■^'3-0.16 

8.86 

6.09 

6.2 

7.6 

5.7 

W^B{Be 

eUe) 

1 yc+O.OS 
' ^-0.04 

2.06 

1.41 

1.47 

1.7 

1.5 

W^B{Be - 


o 00+1.65 

O.^0-1 33 

- 

2.2(5) 

4.71 

- 

4.78 


V. SUMMARY 

In this work, we have analyzed the radiative leptonic Be —)■ jiiy decays in the NRQCD 
effective field theory. NRQCD factorization ensures the separation of short-distance and 
long-distance effects of Be —)■ 'yi9 into all order of ag. Treating the photon as a collinear 
object whose interactions with the heavy quarks can be integrated out, we arrive at a 
factorization formula for the decay amplitude. 

We have calculated not only the short-distance coefficients at leading order and next-to- 
leading order in ag, but also the nonrelativistic corrections at the order |vp in our analysis. 
We found that the QCD corrections can sizably decrease the branching ratio, which has 
very important impact on extracting the long-distance operator matrix elements of Be- For 
phenomenological applications, we have estimated the long-distance matrix elements, which 
are further used to explore the photon energy, lepton energy and lepton-neutrino invariant 
mass distribution. These results can be examined at the LHCb experiment. 
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Appendix A: Ward identities for matrix elements 


In this section, we will derive the constraints on the t 7 form factors following a 
Ward identity for the conservation of the electromagnetic current. To be more specific, let 
us consider the following matrix element: 

( 7 (A;,e)I( 07 ^ 756 )(0)|5c) = iee*^^ J ( 07 ^ 756 )(0)|5c). (Al) 


In this case, the electromagnetic current includes contributions from heavy quarks j®'™' = 
eccy^c + 

The conservation of the electromagnetic current implies a Ward identity for the matrix 
element of the time-ordered product in (Al) 

[ dV-"(0|TA-(x)|(c7,756)(0)|:Be) 


Qn X I x'^^oo 

r0_ 


I —00 


= j dW'''^((0|j^-“-(x) (07^756)(O)| 5 c) 0 (a:°) + ( 0 |(c 7 ^ 756 )( 0 ) j^-“-(a:)|Se)0(-a:°))[ 

= j d^a;e“*^'^((01(07^756)(0)1 Ee) - ( 0 |(c 7 ^ 756 )( 0 ) j^-“-(x)| 5 c)) 

= j . (ct„756)(0)1|S.) . (A2) 

The commutator on the right-hand side is non-vanishing since the operator c'j^'y^b carries 
an electric charge. It can be evaluated as: 

j d^xe-^^'^{ 0 \[j^-’^ix ), (07^756)(0)]|;B(psJ) 

= [ d^xe-^’^'^{0\[eccl^{x)c^{x) + ebbl^{x)bm{x) ,cl{0){-f^-f^j5)nsbsiO)]\B{pBj) 


= (cc - Cb) ( 01 ( 07 ^ 756 ) ( 0 ) 1 5c( pbJ) 
= i{ec- eb)fBcPB^,u ■ 


(A3) 


The most general parametrization of the matrix element on the left-hand side without 
can be written in terms of five form factors fi{k‘^,PBc ' k) 


I j d‘‘ie“-{0|TjJ"' (a:) (<r,„75f>)(0)|Bc) 


i[fig,,u + f 2 PB,,^^PB,,u + fsk^^ku 
+hkf,PB,,u + Upb,M . 


The Ward identity (A3) implies two constraints on these form factors 

{Pb, ■k)f2 + k^fi = {Cc - eb)fB ,, /i + k'^fs + {Pb, ■k)fb = 0. 


(A4) 

(AS) 


For a real photon k"^ = 0, these constraints fix uniquely the form factor f 2 ( 0 ,pb^ ■ k), and 
relate /i(0,pb^ ■ k) and / 5 ( 0 ,pb^ ■ k), which leads to 


( 7 (e, A:) 1 07 ^ 756 15c( pbJ) = iepB, ■ kf^ j^ fscPB^f^PB, ■ e*. (A 6 ) 

This is the same as the result in Eq. 0 as presented in text, with the identification PBpkf^ = 

A. 
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Appendix B: Passarino-Veltman integrals 


The coefficients c* and di are related to the scalar Passarino-Veltman integrals dehned 
in Ref. [7^1 ETj, and we have split the hnite pieces bi = Ci = and di = 

Bi = Bo (O, z'^ml, z^ml) , 

B 2 = Bo , 

Bo = Bo {ml {y^ - zz) /z, 0 , ml) , 

Ba = Bo{y^ml,ml,z^ml) , 

Cl = Co {ml, 0, ml {y^ - zz) /z,0, ml, ml) , 

C 2 = Co {z‘^ml,y^ml,0, ml, z‘^ml, ml) , 

Co = Co {ml, z‘^ml,z‘^ml, ml, 0, z^ml) , 

Ca = Co {ml {y‘^ - zz) /z,mly‘^,mlz^,Q,ml,mlz^) , 

Di = Do {ml, z‘^ml, y^ml, 0, z'^ml, ml {y^ - zz) jz, ml, 0, z'^ml, ml) . (BI) 


Here we give the the results of divergence integrals. 


Ri = — -fin ^ 


^uv 


9 2 ’ 


euv mt 


euv ml y^ - zz 


Ba = —+ ln^^ + 2 + ^(7i(j/)ln(^**'^^ 
ec/v y^ml ^ 


2 = 1 
,2 


= 

Di = 


1^1 . h' ^ 

-h ti -|- In — K — 2 — 


2zml eiR 
z 


mt 


ii{y) 

2lnz 


) -ln(7*(2/) - 1)), 


2mfz{Pml z {y - z){y + 

-2y^ In y - {y‘^ + z‘^ - l)\n z - y‘^{l + 2\n2)) + {-go + y"^ + z"^ - l)gi 
+ {95 + ~ + 1 ) 5'2 + {—95 + y'^ ~ Z + l)go + {90 + y"^ + z"^ — 1 ) 5 ' 4 )), 


(B2) 
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where 


71,2(3:) 

9i 


92 

93 

94 

95 


(a ;2 — + 1 )^ — 4x^ + + 1 

In (y^ — z'^ + 1)^ — 4y^ — — z^ + , 

In ^-y/ (y^ — z'^ + 1)^ — 4y^ + y^ — z^ + , 

In (y^ — z^ + 1)^ — 4y^ + y"^ — z^ + 1^ , 

In (y^ — z^ + 1)^ — 4y^ — y^ — z^ + 1^ , 

\ly"^ — 2y‘^ {z"^ + 1) + (yz"^ — Vf. 


(B3) 


Appendix C: One loop corrections to the axial-vector form factor A 

The most general structure of the matrix element of the axial-vector current is 
parametrized by: 

( 7 (e, A;)|c 7 ^ 756 |[c 6 (^ 4 ^^)]) = ie (Cl) 

This section will be devoted to demonstrate the gauge invariance at the one-loop level in 

as, namely 

A" = = A, (C2) 

15^ = 15. (C3) 
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The contributions from individual diagrams to are given as 


Af = 


ebCpasy/^c _ 4:{z - l)z'^ y^z"^ - 2{z - l)z^ ~ i 1 /^ ^ 

47 rmb _iiR [y-i - pf {y^ - pf z ^ {y"^ - P)z 

2y^z -y^ + {z^ - 1 )^ -2y^z + {z - 1)P + y"^ 

z{y^-Pf " z{y^-Pf " 

2{‘iz - 1)P + y^-y^ {z^ + 4;^ + 3) ^ {z^ - 1) {y^ - z^ + 1) 

{P-y‘^)z {y^-P)z """ 


+ (— y"^ -\- z^ -\- Az — l) di , 


= 


efeCira^y2]Vc ri 1 , z + y"^ 


^ ^ \ ^ h J- 

Anrrib 2 iuv y'^ — P ^ 


-y^{2z + 3) + 2z^ + 5z^ + 2z-l. 
2z {P — y‘^) 


y^ + z^ + Az + 3 ^ ^-3y^ + z^-l), , 2,2 


y2 _ p 


+ 2z{y^-P) + +^-l)c 4 j, 


— 

— 


A! = 


CfeCira^y^l'l 1 y’^-z^W y^-zz y^-z^\\ ^ ] 

-^^- ~^~2 -^^ 2 + -^ 77 ^- z^h-zci , C 6 

47rmb [2 euv U ~ z^ y^ — z^ 2 [y^ — z^) 


ebCpUsV^c r 1 , - z"^ + 1 


— 2:^ + 1 . 


P -^ oTT-A“ A TT^—— 771 -T'^s , (C7) 

2euv 2 (|/2 — zz) 2y^ — 2zz 2 [y^ — zz) 


The mass counter term and wave function renormalization give the contributions: 

^CT-m = 0) 

^CT-F = — VcT-F- 

The contributions from individual diagrams to A^ are given as 


^ ebCFasV2Nc\ z {y"^ {-1 z"^ + IQz + 1) z^ + {z - 1)P + y"^ {^ly"^ + ?,z^ - ^z - ll)) 


K = 


Anrub 


(y^ - zz) {y^ - ypy 


2y‘^{‘i-2z)P + zP+ y\-iz + 2))^ z + y^ , 2z {-y‘^ + z^ + 3) 

H-77-177--T^^2-FT-7T-C4 

(l/3 - yz^)^ y^ - y^zz [y^ - zy 

—y^P {y^ + z'^ — Az + 3) + {z^ — Sz"^ + 5z + 1) P + y^ 

z {—y^ + z'^ + z) {y"^ — py 

^y‘^{3z-h)P - {z-l)P+ y^{y'^ + z^-l)^] ^ 

2 : (|/3 - ypf \ 

f. ebCFC(sp2NcP {—Tz"^ + lOz + 1) P + (z — 1)P + y^ {3y‘^ + 3z‘^ — 3z — 11)) 


A," = 


47rmfe (y2 _ 2 ;^) _ ypy 

2y‘^{2z-3)P - zP+ y^{z + 2) z + y^ 2z {-y^ + z^ + 3) 
(^3 _ ypy y'^- y^zz ^ (y"^ - P) 

—y‘^P iy"^ + 2 :^ — 4^: + 5) + {z^ — 3z'^ + 5z + 1) P + y^ ^ 
z {zz — y"^) {y"^ — py 

y‘^{3z - 3)P -{z- l)P + y"^ {y"^ + z"^ - 1) ^ ^ ^ ^ ^ 

-—- —2 -O 4 J + 

2 : ( 7 ^ - yz^) 


(CIO) 
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K = K + 


ebCpasV^^rZ {by^ — — 6z — 1) z (3?/^ — 3^^ — Az — 1) 


ATlTTlb 


{y^ - z"^) {y^ - zz) 


+ 


z{-by^ + bz^ + Az + l) 

^ (|/2 - P ) (|/2 - zz ) ^ ’ 


= A! + 


ebCpOLgy/^Nc r bz 1 


+ 


(|/2 - P) (|/2 _ 


+ 


(Cll) 


Airmb y"^ — z^ iuv —y"^ + z"^ + z iy"^ — z^) (y^ — zz) 
z (3|/2 - 3^2 -4z-l)^ 


-bs. 


(y^ - z^) (y^ - zz) 

Similar, the mass counter-terms and wave function renormalization corrections give: 


(C12) 


^CT-m 


ebCpasy/^c r 3g 1 4 

Airrub z^ — y"^ iuv 3 ’ 


— k’^ 

^CT-F — ^CT-F- 


(C13) 


Adding the above contributions, one may derive the relation A*^ = A^, which is guaranteed 
by gauge invariance. One can obtain the one-loop results for A by adding up the anti- 
symmetrical part with Cb —t Cc and nib t-t nic- 

The contributions from individual diagrams to are given as 


Of = 


^ ebCFOis\/‘2Nc. 2 y'^ {bz^ — Qz — 1) z^ — {z — l)z^ + y'^ {—‘iy‘^ + z"^ +'8:Z + 7) 

An iiF 1/2 (^y2 _ ^y2 _ 

' o~ / o ~o\ ^1 ~r o~ / o ~\ ^2 


y^z (|/2 — p) 


y‘^z (//2 — zz) 


Ay'^z + {z'^ - 4 ^: - 1 ) ^2 _ ^4 _ I'^p + y'^ + 2 y‘^ {z‘^ - z - 2 )) 

-03 H-o / o- bi 


Z (|/2 — P) (//2 _ zz) 
-2 zc4 -h 4 ^: 03 ], 


y2z (jj2 _ p^ 


{CIA) 


. ebCpasV^Jfcr 1 y"^ (-5^2 4z + l)z^ + {z - l)z^ -f {y"^ + bz"^ -2^-5) 


15^ = 


\.~l ^ 

^uv 


An 

z {zz^ + y ‘^{2 - 3z)) 
y^ — y'^P 


y2 (^y2 _ ^ 2 ) (^y2 _ ^ 


, ^2 -Ay'^z -F (- 2:2 4 ^ _)_ X) ^2 ^ ^4 

bi + - - —b2 + --&3 


1 /^ — y'^zz 


Z (//2 — p) (|/2 _ 2 ;^) 


, z{{z-l)z^-yHy^ + 2z^-2z-A))^ , ^ 

T-2- ^ -2^2- O4 + ^2:04 , 

yZz l^y^ _ 

. _ ebOpgsV^c r 1 {z^ + lOz + 1)2"^ + y^ - 2y‘^ {z^ -H 62 : + 5) 

47 r iuv {y‘^ - z"^) {y‘^ - zz) 

z (—5//2 + 52;2 _|_ 6 ;^ 4 _ 1 ) (^z^ 4 - Qz 4 - 1 ) p 4 - — 2y‘^ ( 2;2 - 1 - 4:2 -f 3) 


(CIS) 


H—^^^^^— -b2 - 

(//2 - p) {y2 - zp 


(|/2 - P) (//2 _ zz) 


bs], (C16) 


. ebCpasV^cp'^ - z'^ - 8 z - 7 1 P-y‘^ z {-y^ + z‘^ - 1) 


Of = 


+ 


+ 


47 r 1/2 _ p y 2 _ (^y 2 _ p-j (^y 2 _ 

(2:2 _|_ g^; -|- 1) ^2 _|_ ^4 _ 2y2 ^^2 _j_ 4^ _|_ 3^ 


(//2 - P) (|/2 _ 2:5) 


-^ 3 ], 


(C17) 
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(C18) 


The sum of them is 


O 


A 

b-e+CT 


2 

>2 

2 

— _A^ 

'^CT-F — ~^CT-F- 


SebCpasV^c {{z - 1) ln(z) - 2z + 2/3ti) 

4:71Z 


(C19) 


We can get the one-loop result in Eq. after adding up the symmetrical part with 

Cft -A- tc and mb O me- 
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